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Abstract 

Recent experimental results on B — > pp decays [1-3] indicate that the CP asymmetry S p + p - will 
give an interesting determination of a = arg [— (VtdY^) / (V u dV* b )]. In the limit when the p width is 
neglected, the B — ► tttt isospin analysis can also be applied to B — > pp, once an angular analysis is 
used to separate transversity modes. The present bound on the shift of S p + p - from the true sin 2a 
is already stronger than it is for S w + W -. We point out a subtle violation of the isospin relations 
when the two p mesons are observed with different invariant masses, and how to constrain this 
effect experimentally. 



* Electronic address: falk@jhu.edu 
^Electronic address: zligeti@lbl.gov 
^Electronic address: yosef . nir@weizmann .ac.il 
^Electronic address: quiini@slac.stanford.edu 



1 



I. INTRODUCTION 



Rates and polarization fractions for various B — > pp decays have been recently mea- 
sured [1-3]. First measurements of CP asymmetries in these modes are expected in the 
near future. This note is a brief comment on the application of isospin analysis to these 
modes, similar to that for tttx channels [4] to extract Standard Model parameters, and in 
particular the CKM phase a = 2 = ar g [ — (^tdYtl) I (^udKT&)]> from these measurements. 
It is important to constrain such CKM phases as precisely as possible in many independent 
ways. Inconsistent results from different approaches could be an indicator of new physics as 
various measurements that are related in the Standard Model can be affected differently by 
possible contributions from physics beyond the Standard Model. Here we comment on the 
need to parameterize the data to allow for the impact of possible 1 = 1 contributions that 
can occur if the two p mesons have different masses. 1 

In the standard parameterization for the CKM matrix, the phase dependence of the 
B — > p % pP decay amplitudes can be written as 

A.. — T -p +i ^ 4- P -p~ if3 

A^j P'ij 6 ~\~ P%j 6 , ( 1 ) 

where A^j describe B + and B° decays, Aij describes B~ and B° decays, and /3, 7 (and 
a = n — (3 — 7) are the angles of the unitarity triangle (for their precise definitions, see e.g., 
Ref. [5]). is dominated by the tree diagram, while comes primarily from so-called 
penguin diagrams. An important role in the CP asymmetries in neutral B decays is played 
by the B° — B° mixing amplitude, which has the following CKM phase dependence: 

M12 = |M 12 |e 2i/3 . (2) 

The dominant CP violating effect in the B — > p + p~ decay comes from the interference 
between the B° — B° mixing amplitude and T + _. As can be deduced from Eqs. (1) and (2), 
this effect is sensitive to the phase a = tc — {3 — 7 (or 2 = tt — 4>i — ^3)- 

The time dependent CP asymmetry in B — > p + p~ can be parametrized as follows: 



r(g hys (t) -> p+p-) - r(B° phys (t) -> p+p-) 

T(B° phys (t) - p+p-) + T(B° phys (t) - p+p-) 



= SV_ sin( Am t) - C + _ cos( Am t) . (3) 



1 By p mass we mean throughout this paper the invariant mass of the pion pair from the decay of that p. 
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If |P + _/T + _| were zero, so that a single weak phase dominates the decay and if, in addition, 
the final state were purely CP even, then S+- = sin 2a (and C+_ = 0). A separation of 
final CP eigenstates is possible with angular analysis [6]; as we will see below the data show 
that the decays to charged p's are dominantly longitudinally polarized and thus CP even. 

Other CP violating effects in B — > pp decays arise from the interference between the 
T and P terms in Eq. (1) or from interference between mixing and P amplitudes. These 
effects not only have different weak phase dependences, but also depend on the amplitude 
ratio \P/T\ and the strong phase arg(P/T). These complicate the relationship between the 
measured CP violation and the phase a. For a given final transversity a (see discussion 
below), this more complicated relation can be parametrized as follows: 



In the case of two pions, Gronau and London [4] showed how to use the six flavor-tagged 
B — > 7T7T rates and isospin symmetry to precisely determine a even in the presence of the 
additional CP violating effects. Later work showed how one can use the isospin relations 
to bound the uncertainties in a, even when sufficient data to complete the full analysis is 
not available [7-9]. These methods can be applied also for the decays to two p mesons. 
The current experimental data implies that the B — > pp case will give a better intermediate 
result. 

II. BOSE STATISTICS AND BROAD RESONANCES 

The vector-vector decays of a spin zero B meson can have orbital angular momentum 
L = 0, 1, or 2. Hence, for two vector particles, they include both even and odd CP modes. 
Since the decaying B meson is spin 0, the total spin of the two vector mesons must be equal 
to and oppositely aligned to the orbital angular momentum, L. Thus, in the case of two 
identical vector mesons, such as two equal mass p mesons, independent of the value of L, the 
combined space plus spin wave function of the the two identical vector mesons is symmetric 
under particle exchange. Bose statistics then tells us that, just as in the case of two pions, 
the isospin of the two p meson state must be symmetric under exchange of the particles, 




(4) 
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thereby eliminating any possible 1=1 contributions. 2 

While the above argument is made in terms of the amplitudes of a given L, it applies for 
all L. Thus it is equally valid when applied to the amplitudes expressed in any alternative 
angular decomposition. The set of basis functions for describing the decays used in the 
experimental analyses are labelled by the transversity a — 0, ||,_L of the p mesons (which 
both must be the same since the initial state has spin zero). Thus, from this point on, our 
discussion will be in this basis. Note that once this basis is chosen there is no longer any sense 
in which one can separate the different orbital angular momentum contributions within a 
given transversity-labelled state. Since transversity-labelled amplitudes are a choice of three 
orthogonal angular basis functions for analyzing the decays, they contain the full angular 
momentum information. Thus we have a complete set of amplitudes, A"- = A[B — > (/o*/o , ') £r ], 
where a is the transversity label and i and j are the charges of the two p mesons. The CP 
of a given transversity state is well-defined, in the case at hand the states a = and || are 
CP even, while the a =_L states are CP odd [6]. 

The above arguments for the absence of / = 1 in each transversity state do not apply 
for general four-pion amplitudes. This contribution exists even when two pion pairs have 
the same invariant mass and angular momentum. Indeed the fact that p mesons have 
a significant width reintroduces the possibility of / = 1 contributions even for a pair of 
longitudinally polarized p particles. In each B — > pp event the invariant mass of each p is 
measured, and the two values can differ by an amount of order of T p , or rather by the width 
of the region allowed by experimental cuts on the data. The B — > pp amplitude for two p 
mesons with charges q±, g 2 , masses mi, m 2 and helicities Ai = A 2 can have a part which is 
antisymmetric under the interchange of the values of mi and m 2 , and thus, by Bose statistics, 
this amplitude is also antisymmetric in the combined (space, spin, isospin) wave function, 
thus allowing odd isospin, despite the fact that L = S. In contrast, the dependence of the 
even- isospin amplitudes on the p masses is symmetric under interchange of m\ and m 2 . The 
different isospin amplitudes do not interfere. Our main point in this note is that the fits to 
data should explicitly include the possibility of the odd-isospin contribution in B — > pp. 

The size of the 1 = 1 contribution is a dynamical question; we make no prediction. We 

2 In Section 6.1.2.2 of the Babar Physics Book [5] this argument is given correctly for the L = 0, 2 case. 
However, an incorrect conclusion that isospin analysis is not possible for the L = 1 component is stated. 
Mea culpa, HQ. 
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cannot rule out the presence of / = 1 contributions of order (T p /m p ) 2 in the total rate. The 
fact that this amplitude must vanish for equal p meson masses gives it a distinct distribution 
as a function of mi and m 2 from the leading even-isospin terms. The leading contribution to 
the rate due to the amplitude antisymmetric in mi and m 2 can be parameterized by adding 
to the fits a term of the form 



mi — m 2 



m p 



2 



B p (mi)B p (m 2 2 ) , (5) 



where B p (s) is the Breit-Wigner. This contribution vanishes where the even-isospin contri- 
bution peaks. The 1=1 contributions in the p + p~ and p^p channels are unrelated, while 
there is no such contribution to p°p°. Note that even-isospin contributions of the same form 
are also possible, e.g., from the cross-term in 



a + b 



(mi - m 2 ) 



m 2 p 



2 



B p {m\)B p {ml) . (6) 



We expect a, b and c to be of the same order, so the even-isospin contribution proportional 
to ab could be comparable to the 1 = 1 component. 

The question is whether the extraction of the leading even-isospin amplitudes [the a 2 term 
in Eq. (6)] is sensitive to possible contributions of the form (5). Independent of whether 
the correction term is dominated by the c 2 term of Eq. (5) or the interference of a and b in 
Eq. (6), the stability of the fit for the a 2 term can be tested. If the addition of terms of the 
form (5) causes the value of the leading term to shift significantly then further tests must 
be made to ensure a stable value for the on-peak amplitudes. If adding such a term does 
not significantly change the result for the leading term, then we can be confident that the 
correct on-peak amplitudes have been measured. 

While the 1 = 1 contribution must be positive, the subleading even-isospin contributions 
may have either sign. Thus, even if a fit to the data finds that contributions to the rate of 
the form in Eq. (5) are small, that could still be due to cancellations. Such a cancellation 
would be accidental in either the p + p~ or the p ± p° channels, and it is unlikely to occur in 
both. Thus, if the fits in both of these modes are insensitive to terms of the form (5), then it 
is probably safe to assume that the 1 = 1 contributions are likewise small. But, as we stress 
above, it is not the size of these terms that really matters here, but rather the stability of 
the fit to the on-peak, equal mass, pp contribution, for which the isospin analysis is to be 
carried out. If the fits are sensitive to terms of the form in Eq. (5), then further analysis, 
and probably significantly more data is needed. 



As an alternative to fitting the data including terms of the form (5), one can eliminate 
effects of any contributions of this form by decreasing the width of the p bands, A, used 
in the fit (or imposing a cut on \m 1 — m 2 |). Once the accepted p band is small enough, 
the result will be stable against further reduction in its width, and also against changes 
to the leading fit parameters when a term of the form (5) is added. At present, BaBar 
uses a band 0.52 GeV < < 1.02 GeV [2, 3] whereas BELLE accepts a narrower range, 
0.65 GeV < < 0.89 GeV [1]. The possible 1=1 contamination in the B — > pp signal 
diminishes for A < Y p at least as (A/m p ) 2 . If the extracted values of the rates are stable for 
different values of A that would indicate that the 1 = 1 contamination is small and we need 
not worry further about these types of terms, whereas results that are sensitive to A would 
indicate that there is a contribution of this type that must be more carefully investigated, 
or excluded by taking a smaller acceptance. 

Clearly both the approach of adding parameters to the fit and the approach of narrowing 
the acceptance have a statistical cost. We are hopeful that, even with the present data set, 
one will be able to see that the impact of possible 1=1 terms is not large. If their effect 
turns out to be important, then more data will be needed to eliminate their impact. 



III. ISOSPIN RELATIONS 



For each transversity, a, the even-isospin amplitudes have relationships similar to that 
for the two-pion amplitudes [4], 



4°" _i_ A<r — A°~ 

rrz -I — "T ^00 ~~ ^+0 i 



1 

+ ^oo = ^-o • (?) 

Each of these equations can be represented as a triangle in the complex plane. Note that 
the triangles corresponding to the different transversity states can be different. 

Tree diagrams contribute to both A J = 1/2 and 3/2 transitions to / = and 1 = 2 
final states, respectively. Since the gluon is isospin singlet, penguin diagrams contribute 
only to AI = 1/2 transitions to / = final states. Since the final p ± p° states have no 
1 = component, A a +{i and A°L are pure tree amplitudes. Therefore, \A+ \ = \A°_ \ and the 
relative phase of these amplitudes is 27 [see Eq. (1)]. The two triangles originating from 
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Eqs. (7) for any given a can thus be superimposed with a common base, A°^ Q , if all the A^ 
amplitudes are multiplied by a factor e 2 * 7 . 

Electroweak penguin amplitudes, unlike gluonic penguins, contribute to both AI = 1/2 
and 3/2 and hence cannot be distinguished from the tree amplitudes by their isospin struc- 
ture. Since electroweak penguins contribute to both and in Eq. (1), one impact of 
such terms would be a possible difference between A +0 and e 2l7 A_ . The size of corrections 
that contribute to \A +0 \ ^ \A-o\ can be constrained by measuring these two rates. The 
average of the BABAR [2] and BELLE [1] results is 

-0.09 ±0.16, 

o.96i8;8* . (8) 

These results are consistent with the isospin relationship A T o = 0, though with current 
precision the test is not particularly stringent. Given this, there is residual uncertainty in 
the extracted value of a that is not constrained by the isospin analysis. While the impact 
of electroweak penguins on the extraction of a from the isospin analysis of B — > nn can be 
estimated to be of order 1.5° [10], estimates of other isospin violating effects employ hadronic 
models and range from negligible to less than 5° — 10° [11]. The impact of these effects is 
expected to be similar in B — > pp. Dedicated analyses are warranted, since both the matrix 
elements of electroweak penguin operators and isospin breaking are different for the pp final 
state. At the present level of accuracy it is reasonable to assume that these uncertainties 
are small compared to those that are bounded by the isospin analysis; we will neglect them 
in what follows. 

Once the branching ratios B[B — > (p l p j ) a ] = \Afj\ 2 are measured, one can construct the 
two triangles and use this construction to measure the relative phase between A^_ and 
e 2n A° h _ [4]. This phase is 2S a defined in Eq. (4). It arises from a combination of relative 
weak and strong phases and the relative magnitudes of the T + _ and P + _ contributions, 
none of which can be reliably calculated. Using the two-triangle construction to determine 
25 a , there is a fourfold ambiguity in the value of this phase, coming from the four possible 
orientations of the two triangles relative to their common base. 

Until the flavor-tagged branching fractions, B[B° — > (p°p°) a ] an d B[B° — > (p°p°) cr ], are 
separately measured, one cannot determine 8 a . However, one can bound it. Among the 



A T0 
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three averaged branching ratios (summed over transversities), 

B + - = \(\A + -\ 2 + \A + -\ 2 ), 

1 

2 
1 

2 



2 

S+u = i(|A +0 | 2 + |A_ | 2 ), 



£ o = ]- (\A 00 \ 2 + \A 00 \ 2 ) , (9) 



the first two have been measured and there is an upper bound on the third. This provides 
an upper bound on £>q for any a. It is significantly smaller than the rate for the dominant 
longitudinal mode in the other channels. This allows us to place a significant bound on So, 
using the construction described above. Explicitly, the bound reads [7, 8] 3 



cos 2* > 1 - ^ + ^-- A 7+ir^ • (10) 



2B° 00 {Bl_-2B% + 2B m f 
B% + 4£°_B° 

This bound can be further strengthened if experiments constrain C°_ [9] and Cq . 



IV. CORRECTIONS PROPORTIONAL TO 1 - / 

For both B — > p + p~ and B — > p ± p° , experiments have determined that the longitudinal 
fraction fo is close to 100% [see Eq. (14)]. Thus, even if the experiments do not distinguish 
the asymmetry in the longitudinal mode alone, one can use the total asymmetry to constrain 
the longitudinal asymmetry. Since we already know from the data that the decay is almost 
purely longitudinal, the correction is small, of 0(1 — f ). Using S+- = J2 a f <?£>+- and C + _ = 
Z)o-/o-C+_, the differences between the transversity-summed CP violating asymmetries and 
those in the longitudinal mode are given by 

= (1 ~ fo) (V ~ g " - * - (/II ~ h) S +~ ~ ^ , 

cj_ - c + _ = (i - /o) - cl - + 2 ci - ) - (/„ - /,) g " - = ^ • (11) 

The S°_ and asymmetries in each of the transversity channels can in principle be 

anywhere from —1 to +1 subject to the constraints (S^„) 2 + (C^L) 2 < I. Thus, the maximal 

3 The bound in [8], quoted in Eq. (10), is the same as the one in Eq. (2.15) of [7] up to terms of 
0[(Boo/B+o) 2 , {B+-/B+0 - 2) 2 ], where [8] is more restrictive. In [8], the weaker bound in Eq. (2.12) 
of [7], cos2(5o > 1 — 2B®q/B\ Qi is referred to as the Grossman-Quinn bound. 
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deviations of the measured asymmetries from those for the longitudinal modes are 

< (l-/ )(l + |Sj_|), 
ICT-CVI < (l-/ )(l + |C°_|). (12) 

In reality, we expect the error in estimating S9_ to be smaller than this upper bound. To 
zeroth order in |P+_/T£_| we have S+_ = — S+_ = S+_. Consequently, we obtain 

- = (1 - f - /|| + f ± ) + 0[{\ - /o) \P + -/T + _\] . (13) 

One further issue that must be considered is the impact of non-resonant contributions to 
B meson decays to four pions, and that of other resonances that yield the same final state in 
this analysis. These could contribute with opposite CP to that of the dominant longitudinal 
mode. Since the angular distribution given by the decay of a spin-1 longitudinally-polarized 
meson is quite restrictive, the contamination due to all such contributions is effectively 
included in the error of 1 — / , the fraction of non- longitudinal contributions. Thus the 
uncertainty due to these contributions is taken into account by allowing for the uncertainties 
in Eq. (12) when determining the CP asymmetries in the longitudinal mode. 



V. NUMERICAL RESULTS 

The experimental values given by BABAR for the three averaged branching ratios defined 
in Eq. (9) are [2, 3] 

= (27™) x 1(T 6 , (/o) + _ = 0.991 ,; 01 ± 0-03 , 

B +0 = (22.5^;I ± 5.8) x 1(T 6 , (/ ) +0 = 0.97+°$ ± 0.04 , 
Boo < 2.1 x 10~ 6 (90% CL), (14) 

while BELLE obtained [1] 

B +0 = (31.7 ±7.lljj;f) x 10~ 6 , (/ ) +0 = 0.948 ± 0.106 ± 0.021 . (15) 

We take £> + _ = £>°_ and £> +0 = £>° , thus introducing errors of order (1 — f ). These 
are much smaller than the present experimental errors on B + ~ and B + o and therefore can 
be neglected. We use the following averages, based on (14) and (15): 

B° + _ = (27 ±9) x 10~ 6 , 

9 



B° +0 = (26 ± 6) x 1(T 6 , 

Ko = (0-61^) x 10~ 6 . (16) 

The value of B 00 is based on scaling the number of signal events given in Ref. [2] and 
conservatively assuming the efficiency for (/ )oo = U which yields the largest rate [12]. 

The first question to be asked is whether the rates in Eq. (16) are consistent with isospin 
symmetry. Note that in the Bq — > limit, we must have = 2£>° . The central values 
in (16) imply a small £>q but £>° _ ~ £>° , thus the consistency with the isospin constraints 
is limited. Indeed, a statistical analysis [13] of the rates in Eq. (16) finds the goodness of 
the fit is only 24%. Since this confidence level is not extremely small, in the following we 
derive limits on 5 assuming that isospin symmetry holds. 4 Using the isospin constraints as 
coded in [13] and the branching ratios in Eq. (16), we obtain the 90% CL bound: 

cos25 >0.83, (17) 

or, equivalently, 

N < 17°. (18) 

Note that even though the statistical significance of B° + _ — 2B° + + 2£>oo 7^ is small, the last 
term in Eq. (10) does play a role. Had we ignored it, we would have obtained cos2<5 > 0.80. 

It is interesting that the small value of B u o/B +u already puts an upper bound on C + _, 
the measure of direct CP violation. For each transversity component, the isospin relations 
imply, for B% /B% < 1/2, 

The 90% CL bound on Bq Q /B° + _ that can be extracted from Eq. (16) yields, to leading order 
in the small quantity 1 — / , 

< 0.53 . (20) 



VI. CONCLUSIONS 



The present measurements of the rates of the various B — > pp decays already yield 
significant limits on the uncertainty in the extraction of a from the CP violating asymmetry 

4 It was pointed out in Rcf . [3] that the small upper bound on Boo I &+o constrains the penguin pollution. 
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in B° and B° decays to p + p~. Given the large branching fractions of these channels, we 
look forward to an asymmetry measurement in the near future which will determine a with 
interesting precision. To ensure the accuracy of the results it is important to include an 
isospin-1 contribution in the fits to data, as in Eq. (5), constrained to vanish when the two p 
mesons have equal masses. We do not expect the impact of this contribution to be large, but 
it could introduce changes of order (T p /m p ) 2 to the best fit parameters. Once this effect is 
constrained experimentally and the CP-violating quantity S + - is measured, B — > pp decays 
promise to provide the best model independent determination of the parameter a for some 
time to come. 
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